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Abstract. We describe 1 an HMC algorithm for dynamical overlap fermions 
which makes use of their good chiral properties. We test the algorithm in 
the Schwinger model. Topological sectors are readily changed even in the 
massless case. 



1. HMC algorithm for overlap fermions for any number of flavors 

Overlap fermions represent a lattice discretization of fermions with the same 
chiral properties as continuum fermions [1]. Properties of overlap fermions 
are reviewed in [2] (see also [3]). In this contribution we would like to de- 
scribe a Hybrid Monte Carlo (HMC) algorithm for the dynamical simulation 
of overlap fermions, which exploits some of their chiral properties. 

We denote by H {n) the hermitian overlap Dirac operator 7sL>(/i) and 
find D\n)D(n) = H*(p). Since [H^(fjL),j 5 ] = [2, 3] one can split H%(p) 
into two parts, each acting in one chirality sector only, H^(fi) = Hg + (fj,) + 
H%_(ji) where, with P± = ±(1 ±75), 

H 2 0± (V) = 1 -^~P± ± ^~P±e(H w )P ± . (1) 



1 Poster presented at the workshop "Lattice fermions and structure of the vacuum", 
October 5-9, 1999, Dubna, Russia. 
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The non-zero eigenvalues of H 2 (fx) are equal in both chirality sectors 
and hence also their contribution to the fermion determinant: 



det(fl£( M )) = det(fl£») > (2) 

The ' indicates that the zero modes have been left out. 

For Nf dynamical flavors the fermion determinant is thus 

[det(DM)] N f = ^[detiD'i^ff = 
H N fM[det(H' 2 (fi))] N f/ 2 = Li NflQl idet(H> 2 ± (ri)] N f . (3) 

We can use this rewriting to get a Hybrid Monte Carlo algorithm for dy- 
namical overlap fermions for any number of flavors. For each flavor we 
introduce one pseudo-fermion of a single chirality: 

det(fl£( M )) = / d4U±e~ s ' ; S p = 4 _1 <f> ± . (4) 

The choice of the chirality is made such as to avoid zero modes: If the gauge 
configuration at the beginning of the trajectory has non-trivial topology, we 
choose the chirality that does not have an exact zero mode of the massless 
overlap Dirac operator. If the topology is trivial, we choose the chirality 
randomly. To take the zero mode contribution into account, we reweight to 
compute observables 

(O) = (» N fMO) ± I (^ IQI )±. (5) 

Having introduced the pseudo- fermions, doing HMC is straightforward. 
We need the contribution from the pseudo-fermions to the force: 

S -§ = =4(1 - mM^§4 ± ; [Hl ± {»T l <P± = x±- (6) 

We use a rational polynomial approximation for e(H w ) written as a sum 
over poles [4, 5]: 

Straightforward algebra then gives (see also [6]) 

+ Se(H w ) + SH W , l t 

X± jr TT X± *~ CQX±-FfrX± + 2^ c k h kXk±-FfrXk± 



^c k x\ ± H w 5 -^H w Xk±- (8) 
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where we introduced 

Xk± = [H 2 w + bl\- 1 X ±. (9) 

The computation of the force requires thus one additional multi-shift "in- 
ner" CG inversion to obtain the Xk±- 

A few remarks are in order: (1) We anticipate that a straightforward 
HMC for dynamical overlap fermions will suffer even more than with stag- 
gered fermions from difficulties in changing topology due to the existence of 
exact zero modes. By working only in one chiral sector, a change of topol- 
ogy is possible, unimpeded by the fermions, as long as the number of zero 
modes changes only in the opposite chirality sector. (2) Accuracy of the 
approximation of e(H w ) can be enforced by projecting out the lowest few 
eigenvectors of H w , and adding their correct contribution exactly [2] . The 
molecular dynamics evolution of the eigenvector projectors P± in Eq. (1) 
can be included using ordinary first order perturbation theory. However, we 
have not included projections in our dynamical fermion code yet. (3) The 
approximation of e(H w ) used in the molecular dynamics steps need not be 
the same as the approximation of e{H w ) for the Metropolis accept/reject 
step. E.g. an approximation, which is smooth around the origin, can be used 
for the HMD part and the more accurate optimal rational approximation 
with projection for the accept/reject step. 

2. Testing in the Schwinger model 

We tested our HMC algorithm in the Nf = 1 and 2 Schwinger model. We 
first look at time histories of the topological charge Q, determined via the 
number of exact zero modes. We see (Fig. 1) that the topological charge 
changes, even in the massless case. 

We compared our HMC results with fiducial results, obtained by a 
brute force approach (exact diagonalization, and then reweighting with the 
fermion determinant of quenched gauge fields). We notice that the accep- 
tance rate does not drop rapidly and the number of CG iterations does not 
diverge as /j, — > 0. 

3. Conclusions 

The non-zero eigenvalues of H^(pi) in each chirality sector contribute identi- 
cally to the overlap fermion determinant. Utilizing this fact, and separating 
the contribution from the fermion zero modes in non-trivial gauge fields, we 
devised an HMC algorithm for any number of flavors of overlap fermions, 
with changes of topology possible even in the massless limit. The trick 
consists in working in the chirality sector without exact zero modes. 
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N,= l, L=8, m 7 l = 3 (13 = 2. 26354) 
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Figure 1. Time history of Q in the 1 flavor model for several quark masses. 

Preliminary tests in the Nf = 1 and 2 Schwinger model show that 
the algorithm works. The topological charge changes. The algorithm works 
even in the massless case. The acceptance rate does not go to zero or the 
CG count to infinity. Further tests on larger systems and in four dimen- 
sions are needed to better judge the usefulness of the algorithm for realistic 
dynamical simulations. 
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85ER250000 and DE-FG05-96ER40979. We would like to thanks the orga- 
nizers for the opportunity to present this poster during the workshop. 
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